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Abstract 

The data of inelastic 16 0+ 16 scattering to the lowest 2 + and 3~ excited states of 
16 have been measured at E lah = 250,350,480,704 and 1120 MeV and analyzed 
consistently in the distorted wave Born approximation (DWBA), using the semi- 
microscopic optical potentials and inelastic form factors given by the folding model, 
to reveal possible refractive structure of the nuclear rainbow that was identified 
earlier in the elastic 16 0+ 16 scattering channel at the same energies. Given the 
known transition strengths of the 2f and states of 16 O well determined from the 
(e, e') data, the DWBA description of the inelastic data over the whole angular range 
was possible only if the absorption in the exit channels is significantly increased 
(especially, for the 16 0+ 16 0"2+ exit channel). Although the refractive pattern of the 
inelastic 16 0+ 16 scattering was found to be less pronounced compared to that 
observed in the elastic scattering channel, a clear remnant of the main rainbow 
maximum could still be seen in the inelastic cross section at E^b = 350 — 704 MeV. 
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1 Introduction 



The study of elastic, refractive scattering of light nucleus-nucleus systems, like 
16q_|_16q an j 16q ) 12q_|_12q^ ^as proven to be very helpful for our understand- 
ing of the interaction between heavy ions (HI) [1]. The main reason is that the 
refractive structures or nuclear rainbow observed in the nucleus-nucleus elastic 
scattering can provide information on the HI optical potential at small inter- 
nuclear distances. Although most of the studies have been concentrated on the 
elastic scattering (see, e.g., Ref. [2] for the 16 0+ 16 system). The refractive 
effects have also been found and discussed for one-neutron transfer reactions 
in the 12 C+ 12 C, 13 C+ 12 C [3,4] and 16 0+ 16 systems [5]. For example, in the 
i6q_|_i6q S y S t eiI1) where the elastic scattering has shown a pronounced nuclear 
rainbow pattern [2], a clear remnant of the rainbow pattern was identified in 
the 0l7ol5 transfer cross section for the ground state transition to 15 1 / 2 - 
[5]. Beside the transfer channel, the inelastic scattering channel of a refrac- 
tive nucleus- nucleus system can also have a well defined rainbow pattern [6]. 
In particular, it has been shown recently by Michel and Ohkubo [7] that the 
transparency seen in some light HI in the form of Airy structure in the elas- 
tic scattering channel could have similar pattern in the inelastic scattering 
channel. 

To complete our study [2,5] of the refractive structure in the quasi-elastic 
i6q_|_i6q sca ttering, we present in this work the results of our analysis in the 
distorted wave Born approximation (DWBA) of the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering data which have been measured in parallel with the elastic 16 0+ 16 
scattering at £ lab = 250,350,480,704 and 1120 MeV [8,9,10,11]. These data 
were taken for the transitions to the lowest 2 + and 3~ states in 16 at 6.92 
MeV and 6.13 MeV excitation energy, respectively, and cover about the same 
wide angular range as that covered by the elastic scattering and one-neutron 
transfer data. Since the transition strengths of the 2f and 37 states of 16 
have been accurately determined from the (e, e') data [12,13], they can serve 
as constraints for the deformation lengths used to generate the inelastic form 
factors for our DWBA analysis. As a consequence, the measured inelastic scat- 
tering data provide a good database for the study of refractive features in the 
inelastic 16 0( 16 0, 16 0') 16 0* scattering. 

A brief survey of the measurements of the inelastic 16 0+ 16 scattering to the 
2+ and 3f states of 16 is presented in Sect. 2. The basic ingredients of the 
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semi-microscopic folding approach to calculate the optical potential and in- 
elastic form factors for the optical model (OM) and DWBA analyses are given 
in Sect. 3. The DWBA results of the analysis of the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering data at £i ab = 250 — 1120 MeV are discussed in Sect. 4. A summary 
and the main conclusions of this work are given in Sect. 5. 



2 Inelastic 16 0( 16 0, 16 0') 16 0* data 

The experimental spectra have been taken at the incident energies E\ ah = 
250 - 480 MeV [8,9,10] using the Q3D magnetic spectrograph at HMI, and 
at 704 and 1120 MeV [10,11] using the SPEG spectrograph at GANIL. The 
spectral lines of the 3~ state at 6.13 MeV and 2 + state at 6.92 MeV were 
observed in these spectra as the strongest inelastic transitions. The method 
used to obtain the experimental cross sections of inelastic scattering from the 
analysis of the inelastic spectra is described in the following. 

Since the projectile and target are identical particles and the same state can be 
excited in the projectile as well as in the target (and also in both by mutual ex- 
citation), every excited state of 16 appears in the spectrum twice at the same 
position of the corresponding excitation energy: (1) a sharp line is observed 
for the target excitation, and (2) a broad line for the projectile excitation (the 
mutual excitation could not be inferred from our data due to a small cross 
section and an increasing three-body background). The line-broadening ob- 
served for the excited states of the projectile is due to the recoil momentum of 
the photon emitted in-flight by the projectile. This effect is directly related to 
the Doppler-broadening usually observed in the 7-ray spectroscopy. The max- 
imum energy shift in the laboratory system is given by AE max = ±(v/c)E~ / 
[14], where E y is the energy of the 7-transition and v is the velocity of the 
photon-emitting outgoing 16 projectile. These maximum shifts can be rather 
large: at the forward scattering angles they are, e.g., for the 7-transition from 
the 2+ state at 6.92 MeV to the ground state A£ max = ±1.45 MeV at the 
beam energy E lah = 350 MeV, and AE max = ±2.45 MeV at £ lab = 1120 MeV 
(angular shifts, which also result from the 7-recoil, are small and not larger 
than the angular resolution). Note that at 1120 MeV, the base width of the 
broadened line is almost 5 MeV wide, and this needs to be compared with the 
experimental energy resolution of about 0.5 MeV at both incident energies. 

In the spectrum of inelastic excitation the lines resulting from the target ex- 
citation are described by Gaussians using the experimental resolution as the 
width, and the projectile excitation by rectangular distributions with base 
widths 2(v/c)E 1 which were folded with the experimental resolution. For a 
given excited state the number of counts for the projectile and target exci- 
tation must be equal. This is required by the symmetry of the system which 
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Fig. 1. Spectrum of the inelastic 16 0( 16 0, 16 0') 16 0* scattering at £ lab = 1120 MeV 
and Gi a b = 3.0° in the excitation energy region of the 3" and 2 + excited states at 
6.13 MeV and 6.92 MeV, respectively. The line shape used for the fit is described 
in the text. 

implies the equal excitation strengths of the same excited state in both the 
projectile and target nuclei at the center-of-mass (cm.) scattering angles O c . m . 
and it — c . m .. Thus, the Gaussian line of the target excitation and the broad 
line of the projectile excitation have been first normalized to the same area 
before they were added. The resulting distributions were then fitted to the 
corresponding excited states in the spectrum by optimizing the normaliza- 
tions. The fit can be performed simultaneously for several different excited 
states of 16 with overlapping widths. Fig. 1 shows a spectrum of the inelas- 
tic scattering to the 3~ and 2 + excited states at 6.13 MeV and 6.92 MeV at 
= 1120 MeV and 6i a b = 3.0°. The number of counts obtained in the fit of 
an excited state represents the sum of inelastic cross sections for the projectile 
and target excitations, which was further divided by a factor of two to obtain 
the single excitation cross section for the DWBA analysis. 

In the analysis of the measured spectra we have used an experimental tech- 
nique, where two-dimensional plots of scattering angle versus position in the 
focal plane are constructed to perform projections on either axis by apply- 
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ing specific gate conditions. To obtain spectra like that shown in Fig. 1 it is 
necessary to set rather narrow angular gates (typically 0.2° at small angles 
and slightly larger at intermediate angles) before the spectra are projected 
to the position axis. This procedure has been applied at the energies of 250, 
350, 704, and 1120 MeV up to the scattering angles of 22°, 16.2°, 7.2°, and 9°, 
respectively, in the laboratory frame. As a result, the angular distributions of 
the 3~ state at 6.13 MeV and 2 + state at 6.92 MeV have been well separated 
in these angular ranges. At larger angles it was not possible to separate the 3~ 
and 2 + excited states using this method. Their sum could be, however, easily 
projected when the gate in excitation energy was set in such a way that both 
broad lines were included within the acceptable energy limits (in the energy 
range where other excited states of 16 were not observed). In this way, the 
angular distribution of the total (2 + + 3~) cross section has been projected for 
these states up to 42.5°, 37°, and 19.1° in the laboratory frame at the energies 
of 250, 350, and 704 MeV, respectively. At 1120 MeV, the range could not 
be extended due to a steep decrease of the cross sections. At 480 MeV, only 
the angular distribution of the sum of the (2 + + 3 _ ) cross sections has been 
projected. 

We also note that the observed spectrum represents overwhelmingly the inelas- 
tic scattering to the 3 _ and 2 + excited states at 6.13 MeV and 6.92 MeV, and 
the strength mixing by the + and 1~ states at 6.05 MeV and 7.12 MeV can be 
neglected. While no peak has been observed for the + state in the angular bin 
covered by our measurement, a slight rise caused by the 1~ state at 7.12 MeV 
could be seen. The dipole excitation strength of the 1" state is however too 
weak compared to those of the 3~ and 2 + states: B(E1 |) 3.4 x 10~ 4 W.u. 
compared to B(E3 |) « 13 W.u. and B{E2 [) « 3.5 W.u. [15], to make a 
sizable contribution to the total inelastic 16 0( 16 0, 16 ) 16 0* scattering cross 
section. 



3 Main features of the Folding + DWBA analysis 

The most important ingredients of the DWBA analysis of the inelastic scatter- 
ing 16 0( 16 0, 16 ) 16 0* are the complex optical potentials (OP) in the entrance 
and exit channels 

[/£> =< 16 O, 16 0|Vf 6 0, 16 O >, E/£J =< 16 O*, 16 0|Vf 6 0*, 16 O >, (1) 

and the transition potential or inelastic form factor (FF) 

U£L =< 16 O*, 16 0|Vf 6 0, 16 O >, with A = 2, 3. (2) 
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For the elastic and inelastic HI scattering, there are two main prescriptions of 
constructing the nuclear matrix elements (1) and (2): 

(I) The phenomenological approach of parameterizing the OP (1) in terms of 
a conventional Woods-Saxon potential whose parameters to be adjusted to 
the best fit of elastic scattering data. The inelastic FF (2) is then given by 
the radial derivative of the OP, scaled by a factor known as the deformation 
length. This prescription for the inelastic FF is known as the deformed-optical- 
potential (DOP) model where the deformation length is normally obtained by 
matching the calculated cross section to the observed inelastic data [16]. 

(II) The (semi-microscopic) folding approach [17] where the real OP and in- 
elastic FF are calculated microscopically using an appropriate effective in- 
medium nucleon-nucleon (NN) interaction and the ground-state and transition 
nuclear densities for the two colliding nuclei. 

While the shape of the inelastic folded potential has a strong dependence on 
the multipolarity A of the transition, the inelastic FF given by the simple 
DOP model (I) has a A-independent shape [18,19]. As a result, the nuclear 
deformation lengths extracted by the DWBA analysis using the DOP model 
can be significantly underestimated [18], especially for high multipoles A > 3. 

In the present work, following the success of the folding potential in the OM 
description of the elastic 16 0+ 16 scattering [2] and DWBA analysis of the 
0l7ol5 one-neutron transfer data at the same energies [5], we use a recent 
version of the double-folding model [20] to consistently calculate the real OP 
and inelastic FF for the DWBA analysis of the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering. Within this double-folding approach [20] the projectile-target in- 
teraction potential in the elastic or inelastic scattering channel is evaluated as 
an energy dependent Hartree-Fock-type potential of the dinuclear system 

Vf= [< *YMu > + < i'j'\v EX \ji >] = V B + V EX , (3) 

i,i'€A 1 ;j,j'eA 2 

where i>d and v ex are the direct and exchange parts of the effective NN inter- 
action. The calculation of Vd(ex) is done iteratively based on a density-matrix 
expansion method [20,21]. In this calculation, the inputs for mass numbers 
and incident energies were taken as given by the relativistically corrected kine- 
matics [22]. Like in our earlier works [2,5], we have used for the effective NN 
interaction ^d(ex) the CDM3Y6 version [23] of the energy- and density depen- 
dent M3Y-Paris interaction. Since the density-dependent parameters of the 
CDM3Y6 interaction were adjusted to reproduce correctly saturation prop- 
erties of the infinite nuclear matter, to use this interaction for a system of 
two finite nuclei one needs to take into account the kinematical transforma- 
tion from the NN frame to the nucleus-nucleus frame. Thus, the CDM3Y6 
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interaction was further scaled by a kinematics modification factor deduced 
from Eq. (19) of Ref. [24]. At the incident energies considered here, this factor 
leads only to a marginal change of the potential strength (from around 1% at 
250 MeV to 5% at 1120 MeV) compared to that given by the same folding 
approach (3) but neglecting the kinematics modification factor. 

To calculate consistently both the optical potential and inelastic form factor 
one needs to take into account explicitly the multipole decomposition of the 
nuclear density distribution p(r) that enters the folding calculation (see details 
in Ref. [20]) 

PjM^j>M'(r) = ]T < JMXfi\J'M' > C xPx (r)[t x Y Xfl (r)]*, (4) 



where JM and J'M' are the nuclear spins and their projections in the initial 
and final states, respectively, C = \/Att and C\—\ for A ^ 0; Px(r) is the 
nuclear transition density of the 2 A -pole excitation. In the present work, we 
adopt the collective-model Bohr-Mottelson prescription [25] to construct the 
nuclear transition densities for the 2 + and 3~ excitations of 16 as 

Px(r) = Sx^P-, with A = 2, 3. (5) 



Here po( r ) is the total ground state (g.s.) density and 5\ is the deformation 
length of the 2 + or 3~ excitation of 16 0. The g.s. density of 16 was taken as 
a Fermi distribution with parameters [26] chosen to reproduce the shell-model 
density for 16 0, and the neutron and proton parts of the 16 density were 
assumed to have the same shape. Such an isoscalar assumption for the g.s. 
density implies [20] that the same nuclear deformation length is used for neu- 
tron and proton parts of the transition density (5). These deformation lengths 
were determined (see Eqs. (3.33)-(3.35) in Ref. [20]) from the measured electric 
transition strengths of the 2f and 3f states of 16 0, B(E2 f) = 40.6 ±3.8 e 2 fm 4 
[12] and B(E3 |) = 1480 ± 50 e 2 fm 6 [13], so that 5 2 = 1.038 ± 0.048 and 
5 3 = 1.825 ± 0.031 fm were used to generate the transition densities (5) for 
the 2f and 3{~ states. Given these inputs for the nuclear densities, the (energy 
dependent) real folded optical potential Vp ^ and transition form factors Vp 2 -* 
and 1/ F (3) were obtained straightforwardly in the double- folding approach [20]. 
Along with the nuclear FF, the Coulomb inelastic form factor Vq^ was cal- 
culated by the same folding method (3), using the g.s. and transition charge 
densities of the two colliding 16 O nuclei and the Coulomb interaction between 
the two protons. 

To have the complex energy dependent optical potentials for the entrance 
and exit channels, the real (folded) OP was further added, as in our previous 
folding analyses, by a Woods-Saxon (WS) imaginary potential that contains 
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both the volume and surface terms. Thus, the OP used in the present DWBA 
analysis is 

u£l ut) (R) = V1 0) (R) + Nfv^(R) - i[W Y (R) + W D (R)}. (6) 



Here is the elastic Coulomb potential between a point charge and a uni- 

1 /3 1 /3 

form charge distribution of radius Rq = 1-3 x (A/ + A 2 ) fm and the volume 
and surface parts of the imaginary WS potential are determined as 

W V (R) = W V {1 + exp[(r - R v )/a v }}-\ (7) 

W D (R)=AW D exp[(R-R D )/a D ]{l + exp[(R-R D )/a D ]}- 2 . (8) 

For uj® , the renormalization factor Ny^ of the real folded potential Vp ^ and 
parameters of the WS imaginary potential were adjusted in each case to the 
best fit of the measured elastic 16 0+ 16 scattering data [2]. It is commonly 
assumed in the DWBA analysis of inelastic nucleus- nucleus scattering [16], es- 
pecially when one uses the DOP model for the inelastic FF, that C/^t = U in ■ 
In the present work, given the highly accurate inelastic 16 0( 16 0, 16 ) 16 0* scat- 
tering data and inelastic folded FF's fixed by the known transition strengths 
of the 2+ and 3j" states of 16 O, we find that this assumption is not always a 
good approximation, and it could lead to an underestimation of the transition 
strength of the nuclear excitation when the refractive effect is strong. 



Within the folding model [20], the following complex energy dependent inelas- 
tic FF is often used for the transition potential (2) 

U£L(R) = vi X) (R) + N^Vj X \R) - ^A^g®, where A = 2, 3, (9) 



and W in {R) = W Y {R) + W D {R) is the WS imaginary part of u£\R). With the 
deformation length fixed by the known transition strength, the renormalization 
factor iVy^ of the real folded FF is usually adjusted to the best DWBA fit 
of inelastic scattering data and it should be close to unity for the folding 
procedure (9) to be meaningful. All the OM and DWBA calculations were 
made using the code ECIS97 written by Raynal [27]. 



4 Results of the DWBA analysis and discussion 



To have a reliable DWBA description of the inelastic 16 0( 16 0, 16 0') 16 0* scat- 
tering, one needs first to determine the parameters of the OP (1) for the 



8 



Table 1 

OP parameters (6)- (8) for the entrance and exit channels used in the DWBA analysis 
of the inelastic 16 0( 16 0, 16 0') 16 0* scattering at incident energies of 250 to 1120 
MeV. Jy and Jw are the volume integrals (per interacting nucleon pair) of the 
real folded potential and of the imaginary WS potential, respectively. If blank, the 
values are the same as those in the preceding row. 



-*— 'iao 


channel 


V 




Wv 




Oy 


Wr, 

' ' lJ 


1 { ] \ 




' VV 


(MeV) 






(MeVfm 3 ) 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


(MeVfm 3 ) 


250 


entrance 
exit (2+) 
exit (3") 


0.781 


290.0 


38.34 
48.46 
34.87 


4.405 


0.949 


6.630 


5.657 


0.483 


98.81 
119.4 
91.76 


350 


entrance 
exit (2+) 
exit (3 _ ) 


0.845 


298.2 


25.32 
35.76 
27.00 


5.832 


0.629 


6.932 


4.846 


0.350 


103.0 
140.8 
109.1 


480 


entrance 
exit (2+) 
exit (3") 


0.791 


261.5 


31.59 
41.89 
33.89 


5.421 


0.711 


1.647 
15.75 
6.586 


5.114 


0.415 


99.88 
161.9 
117.6 


704 


entrance 
exit (2+) 
exit (3 _ ) 


0.867 


256.1 


28.62 
39.96 
32.14 


5.669 


0.652 


7.095 
6.010 
6.862 


4.153 


0.450 


107.7 
144.2 
119.1 


1120 


entrance 
exit (2+) 
exit (3 _ ) 


0.905 


218.6 


46.34 
61.98 
50.30 


4.793 


0.809 


0.868 
3.888 
1.178 


5.640 


0.391 


109.1 
152.7 
119.0 



entrance and exit channels. In addition to the scaling the CDM3Y6 interac- 
tion by a kinematics modification factor, a more accurate local approximation 
[21] was also used for the exchange potential Vex compared to the folded po- 
tentials used earlier in Ref. [2]. Therefore, we have made again a detailed OM 
analysis of the elastic 16 0+ 16 scattering data [8,9,10,11] using the complex 
OP (6). The obtained OP parameter sets (see the best-fit OP parameters for 
the 16 0+ 16 entrance channel at different incident energies in Table 1) re- 
produce the same refractive structure in the elastic 16 0+ 16 cross sections as 
that reported in Ref. [2]. Note that the improved folding procedure leads to 
a slightly smaller renormalization factor of the real folded potential at 
all energies, with the largest difference of about 10% found at 250 MeV. The 
corresponding volume integrals Jy are changed only by about 2-5% and they 
indicate practically the same potential family for the real OP as that discussed 
in Ref. [2]. 
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To show more clearly the evolution of the refractive Airy pattern with the 
increasing incident energy, we have plotted in Fig. 2 the measured elastic 
16 0+ 16 data and the calculated cross sections versus the momentum transfer 
q = 2&sin(© c . m ./2), where k is the wave number of the projectile. In such a 
presentation, the location of the first Airy minima (Al) can be traced up 
to the energy of 704 MeV. The rainbow pattern associated with Al [2] is 
somewhat obscured at = 250 MeV by the Mott interference caused by 
boson symmetry between the two identical 16 nuclei, but it becomes more 
prominent at E^b = 350 MeV. Actually, the elastic data at 350 MeV [9] is now 
known as a very clear experimental evidence of the nuclear rainbow scattering 
pattern in elastic heavy-ion scattering. In the present work we investigate 
whether this rainbow pattern also exists in the inelastic 16 0( 16 0, 16 ) 16 0* 
scattering channels. 

As discussed earlier in Ref. [2], a very consistent absolute angular calibration 
was obtained for the elastic 16 0+ 16 scattering data at E^ = 1120 MeV, 
where all the gauging conditions were well fulfilled for the target constituents 
6 Li, 12 C, 16 0, 40 Ca and 51 V simultaneously. The inelastic 16 0( 16 0, 16 0') 16 0* 
scattering data were also quite accurately measured at this energy, with the 
inelastic scattering cross sections for the 2 + and 3~ states of 16 well sep- 
arated over the whole angular range covered by the elastic scattering data. 
Since the DWBA formalism also works better at medium and high energies 
(where the coupled-channel effects are not significant [16]), we started first our 
DWBA analysis of the inelastic scattering data at 1120 MeV, and the DWBA 
results are compared with the measured 16 0( 16 0, 16 0') 16 0* angular distribu- 
tions in Fig. 3. Given the deformation lengths fixed above by the measured 
electric transition rates of the 2 + and 3~ states and the validity of DWBA, the 
complex FF's (9) of the inelastic scattering to these states were expected to 
give some reasonable description of the data with the unrenormalized folded 
FF (iV v 2) = = 1). However, the DWBA calculation based on such an 
assumption, using the same complex OP in the entrance and exit channels, 
turned out to grossly overestimate the measured angular distributions for the 
2 + and 3~ states of 16 (see dotted curves in Fig. 3). Even if we assume lower 
values for the deformation lengths, 5 2 (3) $2(3) — ^2(3) where A<5 2 ( 3 ) are the 
experimental uncertainties of the adopted 5 2 (3) values discussed in Sec. 3, the 
absolute 2 + and 3~ cross sections are reduced only by around 10% and 4%, re- 
spectively. Such a small change in the inelastic cross section is almost invisible 
in the logarithmic scale of the plot. 

We tried further to adjust the renormalization factors of the real folded FF's to 
best fit the data, using the same imaginary inelastic FF's. These factors then 

(2) (3) 

turned out to be significantly smaller than unity: N y ps 0.54 and N y m 0.06 
which are too unrealistic. Moreover, with iVy approaching almost zero, the 
strength of the imaginary 3~ inelastic FF is still so strong that the DWBA 
cross section is roughly twice the measured cross section over the whole angular 
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range. Keeping Ny^ value within the range of 0.7-1.0, which is close to Ny 3 ^ 
values obtained in numerous folding analyses of elastic a-nucleus and nucleus- 
nucleus scattering, a reasonable description of the data could only be reached 
by switching off the imaginary part of U^l ns (see dashed and dash-dotted 
curves in the lowest part of Fig. 3 which were obtained with Ny^ = 0.938 and 
1, respectively). For the 2 + excitation, even after the imaginary part of U^l ns is 
neglected, we could not get a good description of the measured cross section for 
the 2 + state over the whole angular range. If the real folded FF is renormalized 
by Ny ~ 0.63 (a strong quenching of the transition strength), the DWBA 
result fits reasonably the data at large angles but underestimates significantly 
the 2 + data at forward angles (dashed curve in the middle part of Fig. 3). When 
the real folded FF is kept unchanged (Ny = 1), the DWBA calculation can 
only fit the data points at smallest angles and strongly overestimate the data 
at large angles (dash-dotted curve in the middle part of Fig. 3). An overall 
good description of the measured 2 + cross section could finally be reached 
only when a more absorptive OP is used for the 16 0+ 16 2 + exit channel, 
with the absorption strength around 40 to 60% stronger than that in the 
entrance channel (see Table 1). In this case, we have used the real folded 
FF with A^ 2,3) = 1 and the same WS geometry of the 16 0+ 16 0* imaginary 
OP as that in the 16 0+ 16 entrance channel, only the WS strengths of the 
imaginary OP in the exit channel were adjusted to obtain a good description of 
the inelastic scattering data by the DWBA calculation. For the 3~ excitation, 
such a procedure lead only to a slightly more absorptive OP in the 16 0+ 16 3 - 
exit. The situation remains nearly the same when we switch on the 3-channel 
coupling 2 + <-> + <-> 3 _ . The coupled-channel (CC) calculation using the 
same OP and inelastic FF as those used in the DWBA calculation leads to 
a slight reduction (around 10-30%) of the inelastic cross sections at largest 
angles (compare solid and dash-dotted curves in Fig. 4). After the depths 
of the WS imaginary potentials were re-adjusted to fit the data in the CC 
calculation (the dotted curves in Fig. 4), the absorptive strength of the OP is 
reduced by only around 5-6%. Such a weak CC effect shown in Fig. 4 confirms 
the conclusion made based on the DWBA results plotted in Fig. 3. 

We have further analyzed the 16 0( 16 0, 16 0') 16 0* data at the lower energies. 
Since the excitation energies of the 2 + and 3~ states of 16 are separated 
by about 0.8 MeV, the inelastic cross sections for these two states could not 
be well separated at large angles due to stronger kinematical effects. There 
only the sum of the 2 + and 3~ cross sections could be measured. We had 
to maintain, therefore, in each case a good DWBA description of both the 
2 + and 3~ data at forward angles as well as their sum at large angles, in 
the refractive angular region. Our DWBA analysis of the considered inelastic 
16 0( 16 0, 16 0') 16 0* data at lower energies shows consistently that the inelastic 
scattering to the lowest 2 + and 3 _ states of 16 O exhausts mainly the strength 
of the real part of the inelastic FF (9). 
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If we consider the imaginary part of the inelastic FF as originating from a 
'dynamic polarization' of the transition potential, similar to that of the mi- 
croscopic OP according to Feshbach [28], then these results would indicate a 
dominance of the direct (one-step) inelastic scattering process and the contri- 
bution from higher order terms to the 16 0( 16 0, 16 ) 16 0* reaction is negligible. 
A similar conclusion has also been made recently by Katsuma et al. [29] based 
on their detailed CC analysis of the elastic and inelastic 16 0+ 16 scattering 
at 350 MeV. 

The DWBA results obtained with the inelastic FF consisting only of the real 
folded FF are shown in Fig. 5. Since a renormalization of the real folded FF 
usually could not deliver a good DWBA description of the data over the whole 
angular range and given the dominance of the direct one-step scattering pro- 
cess, we have used throughout our DWBA analysis Ny^ = 1 and 5\ values 
given in Sec. 3. These DWBA results confirm again that a significantly more 
absorptive OP must be used for the exit channels, especially for 16 0+ 16 02+ 
where the absorption was found roughly 40-50% stronger than that in the 
i6q_|_i6q en t rance channel (see the volume integral Jw of the absorptive WS 
potential in Table 1). We note that the enhancement of the absorption found 
for the 16 0+ 16 0* exit channels is very similar to the enhanced absorption 
established earlier in our (finite-range) DWBA analysis of the one-neutron 
transfer 0l7ol5 reaction [5], where the absorption in the 17 0+ 15 exit chan- 
nels is also about 50% or more stronger than that in the 16 0+ 16 entrance, 
with the corresponding volume integrals Jw of the absorptive WS potential 
close to those normally observed for a HI system with strong absorption. 
Since the enhanced absorption implies a shorter mean free path [30,31], our 
DWBA results for the inelastic 16 0( 16 0, 16 ) 16 0* scattering also support the 
conclusion of Ref. [5] that the nuclear mean free path is quite different as 
one considers different reaction channels. The nuclear mean free path in the 
overlapping region was shown [5] to gradually decrease when one goes from a 
tightly bound and spherical double closed-shell nucleus 16 to a less bound 15 
nucleus (in its ground and excited P3/2 states). Now we see the same trend for 
the 16 nucleus in the excited 2 + and 3~ states. Thus, the enhanced absorp- 
tion found in the exit channels should, in general, imply a shorter mean free 
path of the excited (and less-bound) nuclear cluster in the nuclear medium. 
To illustrate the trend of the increased absorption in the exit channels of the 
quasi-elastic 16 0+ 16 scattering, we have plotted in Fig. 6 the volume inte- 
gral of the absorptive WS potential Jw per interacting nucleon pair. While a 
volume integral Jw ~ 100 MeV fm 3 was found for the elastic 16 0+ 16 scat- 
tering channel, with the weakening or even disappearance of refractive effects 
in the inelastic 16 0( 16 0, 16 0') 16 0* scattering and one-neutron 0l7ol5 trans- 
fer channels [5] the correct description of the data over a large angular range 
required consistently an increased absorption in the exit channels. For the in- 
elastic scattering to the excited 2 + state of 16 O and one- neutron transfer to 
the excited 3/2~ state of 15 0, the enhanced absorption has lead to a volume 
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integral Jw up to around 150-200 MeV fm 3 (see Fig. 6) which is close to that 
observed for strong-absorbing HI systems. For example, the OM analysis of 
the elastic scattering of 390 MeV 20 Ne on 12 C using the real WS and folded 
potentials has given J w = 207 MeV fm 3 and 155 MeV fm 3 , respectively [8]. 

The enhanced absorption found for the exit channels of the inelastic 16 0+ 16 
scattering and 0l7ol5 transfer reaction stresses the need to have a realistic 
choice for the OP not only in the entrance but also in the exit channel. The 
use of the same complex OP in both the entrance and exit channels might 
lead to a large uncertainty in the deduced transition strength if one follows 
the standard method of scaling the inelastic FF to match the DWBA results 
to the measured angular distributions. In particular, one needs to be aware of 
this effect while analyzing the inelastic scattering data measured at 'refractive' 
energies for a light HI system containing an unstable nucleus with the excited 
state being either unbound or very weakly bound. We note finally that the 
subtle effect of absorption enhancement could have been found only owing 
to very accurately measured inelastic 16 0( 16 0, 16 ) 16 0* scattering and one- 
neutron 0l7ol5 transfer data which cover a wide angular range and about 6 
orders of the cross-section magnitude. 

We consider now the refractive features in the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering in more details, and the best data set for this purpose are those 
measured at Ei^ = 350 MeV. At this energy, a broad shoulder-like primary 
rainbow pattern preceded by the first Airy minimum at O c . m . ~ 44° (or q m 6.1 
fm -1 ) has been observed in the elastic 16 0+ 16 scattering cross section (see 
Fig. 2). From the DWBA results given by the same OP in the entrance and 
exit channel (dotted and dash-dotted curves) for the inelastic scattering at 350 
MeV shown in Fig. 7 one can see a weak rise in the sum of the 2 + and 3~ cross 
sections at O c . m . > 46°. This broad pattern at large angles is about the same 
for the 2 + and 3~ cross sections (with overlayed oscillatory structures at the 
forward angles which are out of phase) and should be of the refractive nature 
and caused by the same interference mechanism that gives rise to the Airy 
oscillation seen in the elastic 16 0+ 16 scattering [7]. Note that the rainbow 
shoulder is slightly shifted towards large angles because of a small decrease 
of the cm. energy in the exit channel. This structure becomes weaker if one 
uses a more absorptive OP in the exit channels to reproduce the measured 
(2 + + 3~) cross section as discussed above. Nevertheless, a remnant of the 
primary rainbow shoulder can still be seen separately in the calculated 2 + 
and 3~ cross sections as well as in the measured total cross section. From the 
calculated and measured 16 0( 16 0, 16 ) 16 0* cross sections at all energies shown 
in Fig. 5 one can trace this weak rainbow pattern in the inelastic scattering 
cross section from the energy of 350 MeV up to 704 MeV. 

At 250 MeV, the nuclear rainbow pattern of the elastic 16 0+ 16 scattering 
was shown to be shifted to scattering angles larger than 90° [2] and there it 
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is invisible due to the strong Mott interference. That is also the reason why 
one can not see it in the inelastic 16 0( 16 0, 16 ) 16 0* scattering. Interestingly, 
the enhancement of absorption found at higher energies is not very obvious at 
250 MeV where the results given by the same OP for both the entrance and 
exit channels deliver also a fairly good description of the measured data (see 
Fig. 8). We have seen the same trend in our recent DWBA analysis of the one- 
neutron 0l7ol5 reaction [5], where the enhanced absorption in the exit channel 
was also not observed at 250 MeV. It is very likely that the enhancement of 
the absorption in the exit channel is a feature associated with the refractive 
nature of quasi-elastic nucleus-nucleus scattering. Such an effect can be well 
expected from the interpretation of the rainbow scattering as a phenomenon 
associated with a weak absorption which allows a deeper interpenetration of 
the two nuclei, with the refractive pattern determined by the nucleus-nucleus 
potential at small distances [32]. 

Given a weak rainbow pattern seen in the inelastic scattering cross section at 
the energies of 350-704 MeV, one might also expect some sensibility of the 
measured inelastic data to the shape of the nucleus-nucleus form factor at 
small distances. To probe this effect, we have constructed the transition form 
factors (2) using both the folding and DOP methods. To obtain the DOP 
form factors, we have first deformed the (best-fit) elastic folded potential by 
the same deformation lengths 5\ as that used to generate the nuclear transi- 
tion densities (5) for the folding calculation. Although the two choices of the 
FF have more or less the same strength at the surface of the dinuclear system, 
they differ from each other substantially at small radii (compare the solid and 
dotted curves in Fig. 9). This difference shows up clearly in the calculated 
inelastic cross section and one can see from Fig. 10 that a reasonable DWBA 
description of the data by the DOP form factors can only be reached when the 
deformation lengths <5 2j 3 used in the DOP model are substantially reduced, es- 
pecially, for the A = 3 case where <5 3 has to be reduced by a factor around 1.67. 
The corresponding transition rate is then reduced to B(E3 f) ~ 531 e 2 fm 6 
which is nearly three times smaller than the adopted experimental value of 
1480 ±50 e 2 fm 6 [13]. Such a deficiency of the DOP form factor, which leads to 
an artificial "hindrance" of the 3~ excitation strength, has been first pointed 
out by Beene et al. [18] and is now confirmed again by our folding + DWBA 
analysis of the inelastic 16 0( 16 0, 16 ) 16 0* scattering. Moreover, with the in- 
elastic scattering data covering a wide angular range, it can be seen in Fig. 10 
that the DOP method completely fails to give a good description of the data 
points at largest angles which are known to be sensitive to the shape of the FF 
at small distances. Like the rainbow pattern in the elastic 16 0+ 16 scattering 
which can be used to probe the real OP at small internuclear distances [32] , we 
have established that the weak remnant of the primary rainbow observed in 
the inelastic 16 0( 16 0, 16 0') 16 0* scattering at 350 MeV is also helpful in testing 
the inelastic FF at small radii. 
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5 Summary 



Inelastic 16 0( 16 0, 16 0') 16 0* scattering data at £ lab = 250,350,480,704 and 
1120 MeV for the transitions to the lowest 2 + and 3~ states in 16 0, covering 
the same wide angular range as that covered by the elastic 16 0+ 16 scattering 
at these energies [8,9,10,11], have been measured and analyzed within the 
DWBA using the semi-microscopic optical potential and inelastic form factor 
given by the folding model [20] . Although the refractive pattern of the inelastic 
16 0( 16 0, 16 0') 16 0* scattering was found less pronounced compared to that 
observed in the elastic scattering channel, a similar evolution of the primary 
rainbow remnant could still be traced in the inelastic scattering cross section 
up to E lah = 704 MeV. 

Given the strengths of the folded form factors fixed by the deformation lengths 
determined from the experimental B(E2 j) and B(E3 |) data [12,13], a rea- 
sonable DWBA description of the measured 2 + and 3~ angular distributions 
has been obtained only if the absorption in exit channels is significantly in- 
creased, especially for the 16 0+ 16 2 + exit, at the energies where the refrac- 
tive (rainbow) pattern was well observed in the elastic scattering channel. Our 
DWBA analysis also shows consistently that the considered inelastic scatter- 
ing to the lowest 2 + and 3~ states of 16 exhausts mainly the strength of the 
real part of the inelastic form factor (9). In terms of Feshbach's formalism 
[28], these results should indicate the dominance of the direct (one-step) in- 
elastic scattering and negligible contribution from higher order terms to the 
16 0( 16 0, 16 0') 16 0* reaction. 

The enhanced absorption found for the exit channels of the inelastic 16 0+ 16 
scattering and 0l7ol5 transfer reaction [5] indicates that the nuclear mean free 
path is decreasing during the transition from the entrance channel containing 
two tightly bound double closed-shell 16 O nuclei in their ground states to a less 
bound exit channel containing either 15 O nucleus or 16 O nucleus in the excited 
2 + or 3~ states. This would be of interest to study this effect in the future 
analysis of the quasi-elastic scattering data measured for the nucleus-nucleus 
systems involving unstable nuclei, where the partitions in the entrance and 
exit channels are very differently bound. 

The weak rainbow remnant observed in the inelastic 16 0( 16 0, 16 ) 16 0* scat- 
tering was shown to be quite helpful in testing the inelastic FF at small radii. 
A comparison made between the folding and DOP models for the inelastic 
form factor shows clearly the failure of the DOP approach in the description 
of the inelastic 16 0( 16 0, 16 0') 16 0* scattering data at 350 MeV, which cover 
both the diffract ive and refractive regions in the angular distribution. The use 
of the DOP model also significantly underestimates the deformation lengths 
of the nuclear excitations, especially for A > 3. Therefore, we recommend 
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again the folding model as a more reliable tool in the analysis of the inelastic 
nucleus-nucleus scattering. 
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Fig. 2. The OM description of the elastic 16 0+ 16 scattering data measured at 

Elab = 250 — 1120 MeV [8,9,10,11] given by the semi-microscopic optical potential 

(6). The OP parameters are those determined for the entrance 16 0+ 16 channel 

(see Table 1). The cross sections are given in ratio to the Mott cross section and 

plotted versus the momentum transfer q = 2k sin(0 c . m ./2), where k is the wave 

number of the projectile. Al indicates location of the first Airy minimum. 
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Fig. 3. The calculated DWBA cross sections for the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering at E^b = 1120 MeV in comparison with the measured 2 + and 3~ inelastic 
cross sections and their sum. The dotted curves are results given by the original 
(-ZVy = 1) complex inelastic FF [see Eq. (9)], dash-dotted curves are given by 
the original real folded FF only, dashed curves are given by the real folded FF 
which were renormalized to the best fit the data, and solid curves are obtained with 
the original real folded FF but using a more absorptive OP in the 16 0+ 16 0* exit 
channels (see Table 1). 
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Fig. 4. The cross sections of the inelastic 16 0( 16 0, 16 0') 16 0* scattering at 
^iab = 1120 MeV in comparison with the data. The solid curves were calculated 
in the DWBA (see explanation for solid curves in Fig. 3); dash-dotted curves are 
given by the CC calculation using the same FF and OP parameters as those used 
for solid curves, but taking into account the 3-channel coupling 2 + «-> + «-> 3~; 
dotted curves were obtained in the same CC scheme, but with the WS parameters 
of the imaginary OP slightly re-adjusted to best fit the data. 
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Fig. 5. The DWBA description of the measured total (2 + + 3~) inelastic 
16 0( 16 0, 16 0') 16 0* scattering cross sections at E Uh = 250 - 1120 MeV given by 

(2 3) 

the oris; inal (Ny ' = 1) real folded FF. The dotted curves were obtained with the 
same complex OP for the entrance and exit channels, and the solid curves were 
obtained with a more absorptive OP in the exit channels (see Table 1). The cross 
sections are plotted versus the momentum transfer q = 2k sin(0 c . m ./2), where k is 
the wave number of the projectile. ^ 
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Fig. 6. Energy dependence of the volume integral per interacting nucleon pair Jw 
of the WS imaginary OP used for the elastic 16 0+ 16 channel, 16 0+ 16 2 + and 
16 0+ 16 3 - exit channels of the inelastic 16 0( 16 0, 16 0') 16 0* scattering, and the 
exit channels of the one-neutron 0l7ol5 transfer [5] to the ground state 15 O g . s . and 
the excited state 15 3 / 2 -. The lines are only to guide the eye. 
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Fig. 7. The same as Fig. 3 but for E iab = 350 MeV. 
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Fig. 8. The same as Fig. 3 but for E Uh = 250 MeV. 
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Fig. 9. The form factors for the transition + — » 3 in the inelastic O+ O scat- 

(3) 

tering at 350 MeV. The folded Vp form factor is the solid curve, the DOP curves 
represent the best-fit elastic folded potential deformed with the same deformation 
length #3 as that used to generate Vp (dotted curve) or with 5% adjusted to the 
best DWBA fit to the inelastic scattering data (dashed curve). 
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Fig. 10. The calculated DWBA cross sections for the inelastic 16 0( 16 0, 16 0') 16 0* 
scattering at E^b = 350 MeV in comparison with the measured 2 + and 3~ inelastic 
cross sections and their sum. The solid curves were obtained with the same folded 

(2 3) 

form factors Vp as those discussed in Fig. 7, dotted curves were obtained with 
the DOP form factors built upon the same deformation lengths £2,3 as those used 
to generate Vp , and dashed curves are the DOP results but with #2,3 adjusted to 
the best DWBA fit to the data. 
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